Toroidal, compression, and vortical dipole strengths in 124 Sn 
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The toroidal, compression and vortical dipole strength functions in semi- magic 124 Sn (and partly 
in doubly-magic 100,132 Sn) are analyzed within the random-phase-approximation method with the 
SkT6, SkI3, SLy6, SV-bas, and SkM* Skyrme forces. The isoscalar (T=0), isovector (T=l), and 
electromagnetic ('elm') channels are considered. Both convection j c and magnetization j m nuclear 
currents are taken into account. The calculations basically confirm the previous results obtained for 
208 Pb with the force SLy6. In particular, it is shown that the vortical and toroidal strengths are 
dominated by j c in T=0 channel and by j m in T=l and 'elm' channels. The compression strength is 
always determined by j c . Is also shown that the 'elm' strength (relevant for (e,e') reaction) is very 
similar to T=l one. The toroidal mode resides in the region of the pygmy resonance. So, perhaps, 
this region embraces both irrotational (pygmy) and vortical (toroidal) flows. 



INTRODUCTION 



The toroidal and compression flows represent a sub- 
ject of intense investigation for many years, see e.g. the 
review The toroidal mode (TM) is known to be de- 
termined by the second order correction to the leading 
long- wave part of the transition EX operator 0-S] • The 
compression mode (CM) being different from TM 

by construction, is related, nevertheless, to TM and also 
represents the second-order correction 0, Q . Despite the 
long studies, the vorticity of the modes is still a subject 
of discussions 0]. Besides, a possible coexistence of the 
dipole TM and pygmy mode (PM) becomes actual 0. 

In hydrodynamics (HD), the vorticity is defined by the 



1- 



curl of the velocity field [9| . Then TM is vortical and CM 
is irrotational However, nuclear models deal with 
the nuclear current rather than the velocity field, which 
causes alternative definitions of the vorticity. In fiol ]. 

the nuclear current component jxi=x+ii r ^ ar i sm g m the 
multipolc decomposition of the current transition density 

< f\jnuc(r)\i >, is treated as unconstrained by the con- 
tinuity equation and thus suitable as a measure of the 
vorticity. Following this definition, the TM and CM are 
mixed flows with both vortical and irrotational elements. 
In our recent study [3] , the relevant vortical operator was 
derived and related by a simple manner to the CM and 
TM operators. The vortical, toroidal and compression El 
strengths were numerically explored with the separable 
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Pb 



random-phase-approximation (SRPA) [l£_ _ 
with the force SLy6 [3]. 

In this paper, the toroidal, compression and vortical 
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El strengths are investigated in semi-magic 124 Sn within 
a wide range of Skyrme forces (SkT6 [l5j], SV-bas fl(ij |. 
SkM* [13, SLy6 [hf, SkI3 0). In addition to T=0 and 
T=l channels, the electromagnetic ('elm') strength rele- 
vant for (e,e') reaction is inspected. The trends with the 
neutron number, from 100 Sn to 132 Sn, are considered. A 
possible relation of TM, CM, and PM is discussed. 



II. MODEL AND CALCULATION SCHEME 

The connection of the vortical, toroidal, and compres- 
sion transition operators with the standard electric mul- 
tipole operator M(EX/j,,k) is discussed in detail in Q. 
The toroidal operator M(tor, A/x) appears as a sec- 
ond order term in the decomposition of M(EXfi, k) in 
terms of momentum transfer k: 



M{EXfi; k) = M{EXjj) + k M(tor, A/x) 



where 



M(EX^) = 
M(tor, A/i) 
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are long- wave limits of the standard electric and toroidal 
operators, respectively; Y\i tl (r) is the vector spherical 
harmonic, and p(f) is the density operator. 

Following the concept [lo} . the vortical operator is Q 
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and its transition matrix elements serve as a measure of 
the vorticity for a given excitation. 

The vortical, toroidal, and compression operators are 
simply coupled 



M(vor, A/x) = Af (tor, A/x) + M (com, A/x) , (5) 



where 



M (com, A/x) = — k M(com', A/x) 



2c 2A + 3 



tfV r A+2 y AM (r) V • j nuc (?) (6) 



is the current-dependent compression operator related to 
its familiar density-dependent counterpart Q 



M(com', A/x) = 



1 



2(2A + 3) 



d *rp(r)r^Y^(r) . (7) 



In ©-(0), the symbol j n uc(r) stands for the operator 
of the nuclear current embracing the convectional and 
magnetization parts: 



jnuc(r) = j c (r) + j m (r). 



(8) 



Explicit expressions for j c (r) and j m (r) can be found 
elsewhere, see e.g. [20| . 

In this paper, only dipole transitions g.s. — > I 71 = 1~ 
are considered. Taking into account the center-of-mass 
corrections (c.m.c.) in T=0 channel-the dipole operators 
to be used in the calculations read 
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FIG. 1: (color online) Photoabsorption cross section in 
124 Sn for different Skyrme parametrizations. The SRPA and 
HF+BCS cross sections are compared to the experimental 
data fTl . 



[2J-|26[ giant resonances as well as El strength near the 
particle thresholds [HIH]. 

For all the modes, the strength function reads 



M(vor, 1/x) = - ± ^ / d 3 r J nuc (r) • r 2 1 



S a (El;E) =3^|(^|A/ Q (S10)|*o)| 2 C(^-^) (12) 
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where (r 2 )o = / d 3 r is the ground state squared 

radius. Note that the vortical operator has no the c.m.c. 

The calculations were performed within the SRPA 
[HI, [U[ ■ The model is fully self-consistent since both the 
mean field and residual interaction are derived from the 
same Skyrme functional. Moreover, the residual inter- 
action includes all the functional contributions and the 
Coulomb (direct and exact) terms. The self-consistent 
factorization of the residual interaction dramatically re- 
duces the computational effort while keeping the accu- 
racy of nonscparablc RPA. SRPA has been successfully 
applied to description of electric [H, [2l| - [23j and magnetic 



where £(E — E u ) is a Lorentzian weight with the aver- 
aging parameter A=l MeV; M a (Elfi) is the transition 
operator of the type a = {vor, tor, com, com'}; "Fo is the 
ground state wave function; E v and |"F„) are the energy 
and wave function of the RPA z^-state. 

The T=0, T=l, and 'elm' channels are defined by the 
proper choice of the proton and neutron effective charges 
e° ff p and gyromagnetic factors g° ff p as: 



T = 

T = 1 
el 



e cff = 1 



9n S p=U9n + 9p), (13) 



„cff 



„cff 



-ef = -1, <L = (14) 



p cS — 1 n cS — Cn 



(15) 



Here g n = —3.82 and g p = 5.58 are free neutron and 
proton gyromagnetic ratios; £ m 0.7 is the quenching 
factor. 

For neutrons in the semi-magic nucleus 124 Sn, the zero- 
range pairing forces are used at the BCS level (HF+BCS) 
[29j . More details of the calculations are found in Q . 
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FIG. 2: (color online) Toroidal strength function in 124 Sn for 
T=0 (left panel) and T=l (right panel) channels, computed 
with the total (black bold line) , convection (red thin line), 
and magnetization (blue dash line) nuclear current. Different 
Skyrme forces are used as indicated. 



III. RESULTS AND DISCUSSIONS 

Results of the calculations arc given in Figs. 1-6. 
In Fig. 1, the SRPA accuracy is inspected for the 
photoabsorption in 124 Sn. Unlike the strength func- 
tions (|12[) . the photoabsorption a is computed with 
an energy-dependent averaging parameter A(i?) in the 
Lorentzian weight, for more details see [28j |. A represen- 
tative set of Skyrme forces with various isoscalar effective 
masses is used: SkT6 (m*=l), SV-bas (m*=0.9), SkM* 
(m*=0.79), SLy6 (m*=0.69), and SkI3 (m*=0.58). Both 
SRPA and HF+BCS results are depicted. 

Figure 1 shows a good agreement with the experimen- 
tal data [lj| for all the forces. The best agreement is 
for SLy6. A discrepancy with the experiment around 
the particle emission threshold (9-10 MeV) may be ex- 
plained by neglecting the coupling with complex configu- 
rations, which is expected to be strong in this particular 
energy region [3(1 HH- Note a small energy upshift of 
the double-peak HB+BCS strength with decreasing the 
effective mass to* from SkT6 to SkI3 (this may be ex- 
plained by spreading the mean- field spectra with to*). 
As seen from Fig. 1, the residual interaction drastically 
shifts the strength to a higher energy, neutralizes the 
m* effect, and transforms the double-peak structure to a 
single-peak one. 

Figures 2 and 3 show the toroidal and vortical strength 
functions in 124 Sn for T=0 and 1 channels and the same 
set of Skyrme forces. The cases of the total ©, convec- 
tion, and magnetization nuclear currents are considered. 
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FIG. 3: (color online) The same as in Fig.2 but for the vortical 
strength. 
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FIG. 4: (color online) The compression strength in 124 Sn 
for T=0 (left panel) and T=l (right panel) channels, gen- 
erated by the operator Since j m -contribution is zero, 
only jnuc = jc case is plotted. 



It is seen that, in agreement with the previous SLy6 re- 
sults for 208 Pb [H, the isoscalar (isovector) strengths are 
dominated by the convection j c (magnetization j m ) cur- 
rent. This is explained by the destructive (constructive) 
interference of the neutron and proton j m -contributions 
in T=0 (T=l) channels. 

As mentioned in Sec. 2. the toroidal and vortial 
strengths represent the second-order El corrections to 
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FIG. 5: (color online) Toroidal (left) and vortical (right) 
strengths in 124 Sn, computed with the total (black bold line), 
convection (red thin line), and magnetization (blue dash line) 
nuclear current. The T=0 (upper plot), T=l (middle plot), 
and electromagnetic (bottom plot) channels are considered. 
Only the force SLy6 is used. 



the leading first-order El response (e.g. photoabsorp- 
tion) . These second-order strengths are less collective [4[ 
and, following Figs. 2-3, more sensitive to the Skyrme 
parameterization. It is seen that, though all the Skyrme 
forces give qualitatively similar pictures, the quantitative 
description of the toroidal and vortical strength functions 
(energy centroids, widths, and gross-structure) notice- 
ably depends on the force. In particular, the SRPA en- 
ergy centroids are generally upshifted with m* , like in 
the HF+BCS photoabsorption in Fig. 1. Note that de- 
pendence on the isoscalar effective mass m* takes place 
in both T=0 and 1 channels. 

It is remarkable that the region 6-10 MeV, often related 
to the pygmy resonance, embraces noticeable fractions 
of the toroidal and vortical strengths, especially in T=0 
channel. For some forces, e.g. SLy6, the strengths are 
even peaked there. This hints that the region 6-10 MeV 
may host both the vortical (toroidal) and irrotational 
(pygmy) motions. 

In Fig. 4, the compression strength function in 124 Sn, 
generated by the operator JTTJ), is presented. The CM 
has no contribution from the magnetization current jm 
and is fully determined by the convection current j c [J]. 
It is seen that, like in the previous figures for the toroidal 
and vortical strengths, the CM demonstrates for T=0,1 a 
significant dependence on the Skyrme parameterization. 
The upshift of the strength with m* (though not regular) 
is also visible. 

In Fig. 5, the compositions of the toroidal and vortical 
strengths in T=0 and 1 channels (|T3")) - (|14[) are compared 
to the 'elm' channel (|T5|) relevant for (e,e') reaction. It 
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FIG. 6: The vortical, toroidal, and compression strength 
functions in T=0 channel in 100 > 124 . 132 s n isotopes, calculated 
within SRPA with the force SLy6. 



seen that T=l and 'elm' responses are very similar and 
both mainly determined by the magnetization current. 
Thus we encounter a remarkable case when the electric 
modes are driven not by j c but j m , and this case may in 
principle be observed in (e,e') reaction. 

The isotopic dependence of isoscalar toroidal, vortical, 
and compression strengths is illustrated in Fig. 6. In 
addition to 124 Sn, the neutron-deficit 100 Sn and neutron- 
excess 132 Sn doubly-magic nuclei are considered. The 
neutron excess (skin) over the N=Z nuclear core is zero 
in 100 Sn and essential in 132 Sn. Fig. 6 shows that the 
most strong isotopic effect takes place for CM. There is a 
considerable growth of the CM high-energy branch from 
100 Sn to 132 Sn, which may be related to increasing the 
neutrons participating in the dipole compression. The 
CM low-energy branch located at the pygmy resonance 
region 6-10 MeV grows from 100 Sn to 132 Sn even more. 
This indicates a strong dominance of neutron oscillations 
in the region. Note that the dipole compression and fa- 
miliar pygmy-like (oscillation of the neutron skin against 
the N=Z core) flows are both irrotational in the HD sense 
and their velocity fields look rather similar. So their mix- 
ture in the pygmy resonance region seems natural. Figure 
6 also shows that, unlike the CM, the toroidal and vor- 
tical strengths weakly depend on the neutron skin. Per- 
haps the vortical (toroidal) motion takes place mainly in 
nuclear N=Z core. The remaining minor isotopic effect 
may be related to the coupling between the neutron oscil- 
lations and toroidal internal motion. Anyway, both these 
collective motions share about the same energy region. 
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IV. CONCLUSION 

The toroidal, vortical, and compression dipolc strength 
functions in 10 °. 124 432g n j so topes were analyzed in the 
framework of the self-consistent separable Skyrmc-RPA 
approach [l2|, HH. A representative set of five Skyrmc 
forces with essentially different isoscalar effective mass 
to* was used. The isoscalar (T=0), isovector (T=l), and 
electromagnetic ('elm') channels were inspected. The cal- 
culations generally confirm the previous results for 208 Pb 
obtained for the Skyrme force SLy6. In particular, it was 
corroborated that the toroidal and vortical strengths in 
T=0 (T=l) channels are mainly provided by the convec- 
tion (magnetization) nuclear current. A close similarity 
of T=l and 'elm' channels was established. This means 
the E1(T=1) toroidal mode represents an unusual case of 
the electric collective motion determined by the magneti- 
zation current and this case can in principle be explored 
in (e,e') reaction. 

The comparison of the results for 100 > 124 < 132 Sn shows 



that both low- and high-energy branches of the compres- 
sion mode considerably depend on the neutron excess. 
The enhanced dipole strength in the low-energy region, 
often called pygmy resonance, is strongly correlated with 
the irrotational low-energy compression mode and vorti- 
cal toroidal mode (probably located in the nuclear N=Z 
core). Altogether, the pygmy resonance region shows an 
impressive coexistence of various irrotational and vortical 
flows. 
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